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1 Introduction 



A quantum system at thermal equilibrium can be completely described provided that one 
knows its grand canonical partition function, which can be expressed as a path integral. 
For bosonic systems, an advantageous feature of the path integral approach is that of 
employing commuting functions instead of non-commuting operators. For fermionic sys- 
tems, however, such an advantage is not obvious to hold, as the integration variables are 
also non-commuting. 

In 1980, KuboQ used the path integral approach to calculate the grand canonical 
partition function of the Hubbard model, using the strong coupling limit and performing 
a perturbative expansion in the hopping constant (£). Even though his result is valid 
for any temperature, one does not have the exact coefficient of (3 ( j3 — tj?) of the 
high temperature expansion of the partition function. Since then, improvements on the 
calculation of the high temperature expansion up to order (/3t) 9 for the Hubbard model 
in two and three dimensions have been reported in the literature [Q. 

Recently, Grandati et al. presented a method to calculate the grand canonical par- 
tition function of self-interacting fermions, by writing that function on a lattice and using 
the properties of the Grassmann algebra to calculate its expansion in powers of the cou- 
pling constant. They calculated the first two terms for the bi- dimensional chiral Gross- 
Neveu model, obtaining an analytical result; however, their approach is model-dependent. 
More recently, Creutz[|J used a numerical algorithm to calculate the generating functional 
of a fermionic model, rewritten on a lattice. He applied his algorithm to a unidimensional 
fermionic system involving a thousand grassmannian variables. He pointed out that this 
approach does not have the sign problems that generally hamper the application of the 
Monte Carlo method to fermionic models. 

We do not write the grand canonical partition function of a self-interacting fermionic 
model on a lattice; instead, we present a new method to obtain the coefficients of its 
high temperature expansion in d- dimensions, where d > 1. But for the expansion in 
(3, this method does not involve any other perturbative expansion (such as, say, in the 
couplig constant of the model). In section 2 we present the method, an extension to the 
one used to calculate the grand canonical partition function of the anharmonic fermionic 
oscillator |J , a quantum model in d = space-dimension. In reference [|], the properties 
of the Grassmann algebra were used to calculate the moments of grassmannian gaussian 
integrals. This general result, together with the diagonalization of matrices A acr (cr =f, J,) 
- matrices that appear when the trace of any fermionic operator is expressed in terms of a 
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multivariable Grassmann integral (see eq. (|9]) ) — allows us to develop a general approach 
to obtain analytical expressions for the coefficients of the high temperature expansion 
of the grand canonical partition function of any self-interacting fermioning model in d- 
dimensions, even in the thermodynamical limit. In section 3 we apply the method to 
Hatsugai-Kohmoto model, a simple toy model that was used by Hatsugai and Kohmoto 
to explain the metal-insulator transition. The solution of this model is very simple, 
and does not require all of the features developed in section 2. In section 4 we present 
our conclusions and future applications of the present approach. In appendix A, the 
diagonalization of the matrices A " ", for arbitrary lattice dimension and arbitrary number 
of points in the lattice, is described. 



2 Expansion in the High Temperature Limit and the 
Grassmann Multivariable Integrals 

The grand canonical partition function of any quantum system in the high temperature 
limit can be expanded in terms of (3 as 



Z(/3;n) = Tr(e^ K ) 



oo / i \n 

E — ( r &[K n ] P n , (1) 

n=0 n - 



where K is given by 



K = H — /iN, (2) 

H is the hamiltonian of the system, \x is the chemical potential and N is the total number 
of particles operator. 

The fermionic creation (aj) and destruction (a,) operators can be mapped into gen- 
erators of the Grassmann algebra {fji, rjj} as follows |7], [|, ||: 

d 

af -> fji and a,- -> — , (3) 

drij 

where i,j = 1,2, ••-,jV. The generators of this Grassmann algebra of dimension 2 2A/ ", 
written explicitly as {fji, ■ • • , fj_\f] rji,---, r/^}, satisfy the following anti-commutation re- 
lations: 
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{r)i, rjj} = 0, {rji, rjj} = and {77*, rjj} = 0. 
The trace of any normal-ordered fermionic operator O is[[?j] 



(4) 



Tr[0] = / J] (9®(r7,r7) e*=* , (5) 

J i=i 

where we use the shorthand notation: fj = {fji, ■ • • jfjtf} and 77 = {771, • • • , J7jv}> an d 
(9®(f/, 77) is the kernel of the fermionic operator O in the normal order. (By "normal 
ordered operator" we mean an operator in which all destruction operators are to placed 
to the right of all creation operators.) Naively, it can be said that the grassmannian 
function 0®(fj, n) is obtained by replacing a] — > f/j and aj — > rji in operator O |5|, 0. 

Let us consider from now on the case where the creation and destruction operators are 
caracterized by the indices (£; a), where £ is a <i-dimensional lattice vector (d = 1, 2, 3, • • •) 
and a is the spin component. The components of vector I need not to be orthogonal. This 
lattice vector could equally represent either the space vector x or the momentum vector 
k. If the fermionic operator O is a product of n normal ordered fermionic operators Q, 
we have |J 



r n-l ^ EE ^(^)[^(i>)-r? CT («>+l)] 

Tr[Qi = jnn n%(^M«)e ? -° x 

£ cr=±l a=0 

x Q®(f] a (£ ] 0), n a $ 0)) Q®(fj a (£; l),r, a (£; 1)) x • • ■ x Q®{fj <T (£; n-l), r) a (£; n-l)). 

(6) 

where we define a =|= +l,cr =!= — 1- The Grassmann variables in eq.(^|) satisfy the 
boundary conditions 

rj a (£; n) = — n a {£\ 0) and r] a (£; u) = 0, for v > n, (7) 

with cr = ±1 and £ stands for any vector on the lattice. Eq.(|6|) is still valid for a product 
of n ordered ordered operators, not necessarily equal. 

Relation (|6]) is used to write the terms of the expansion of the grand canonical partition 
function in the high temperature limit as multivariable Grassmann integrals. For a d- 
dimensional fermionic model, the coefficients of the expansion Z(/3,ij,) in eq.(|l]) become 
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n-1 



E E E ^>)MM-« CT (^>+i)] 



Tr[K n ] = / J] II II e^^ 1 ^ 

I* <r=±l a=0 

xJC®{fj a {£; 0),r, a (£; 0)) /C®(^(^ 1), ^1 1)) x • ■ ■ x JC®^ n - 1), r) a {£; n - 1)), 



the boundary conditions (0) still hold for the generators r\ a {£\ v). 

It is much easier to handle generators with one index, then we map the generators 
t] a (£, v) and fj a (£, v) into single-indexed anti-commuting variables. The sum in the expo- 
nential on the r.h.s. of eq.(|8]) can be written as 



n-1 



CT=±1 V=0 



E E E^(Mk(M-^(^+i)] 

'=0 

2nN d 

= E Vi A u Vj- 



(9) 



i,j=i 



Note that the argument of the exponential on the r.h.s. of eq.(||) is diagonal in the indices 
£ and a. Having the components of the column vector rjj (or the line vector fj) grouped 
according to the values of a, and then each subset ordered according to u, and finally 
each subsubset ordered according to £, the matrix A will have the block-structure 



(10) 



' A TT CD 
CD A u 

whose entries are matrices of dimension nN d x nN d and d is the dimension of the vector £. 
The indices /, J are such that /, J = 1, 2, • ■ ■ , 2nN d , where N d is the number of points in 
the lattice. The matrices A^T and A^ are identical. Taking into account the anti-periodic 
condition in temperature (|7|) in eq.(^), the matrices A aa , where a =|, j, are found to have 
the following block-structure : 



A^T = A^ 



I ^.]\fd x ]\fd — ljyd x 7Vd Q)jVd x 7V<i • • • Q)j\fd x ]^d N 

®N d xN d ~^-N d xN d — ~^-N d xN d ' ' ' ®N d xN d 



N d xN d 



'N d xN d 



'N d xN d 



t N d xN d / 
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The symbols t N d xN d and Q^vd xAr d stand for the identity and null matrices of dimension 
N d x N d , respectively. Any lattice vector I can be written as 

£ = £ 1 u 1 +£ 2 u 2 + --- + £ d u d , (12) 

where £i = 1,2, ■ ■ ■ , N (i = 1, 2, • ■ ■ , d), and N is the number of points in the lattice in 
the direction of the <i-dimensional basis vector Uj. A particular basis for which K a<T has 
the block-form shown in (|TT|) , yields the following mapping 

T) c {£;v) > V[(i^l n+u ] N d +ei+ ^ 2 _ 1)N+ ... + ^ d _ 1)N (d-i) ■ (13) 

The generators f] a {£; v) have an analogous mapping. With the newly indexed generators, 
the expression of Tr[K n ] (eq.flD) becomes, 



Tr[K n ] = / H dr li d Vi e 1 '^ 1 x 
"* i=i 

xK®{fj, ri;u = 0) K®{fj, m v = 1) ■ • ■ }C®(fj, rj; v = n - 1). (14) 

Note that expression ([14]) , up to the constant — h, is the coefficient at order [3 n of 
the expansion in the high-temperature limit of the grand canonical partition function for 
any self-interacting fermionic model. The specific model to be studied is represented by 
the grassmannian function /C®, but the matrix A is the same for all fermionic models. 
Once the sub-matrices A^ and A^ are null, the multivariable integral ( [14] ) is equal to 
the product of the contributions coming from the sectors: a a =|f and a a ={{ separately. 
The Grassmann functions JC® are polynomials in the generators of the algebra. Therefore, 
the r.h.s. of eq.(I4~) are moments of the multivariable Grassmann gaussian integrals. In 



reference |§ it is shown that these integrals can be written as co-factors of the matrix A. 
The integrals in eq.(0), for sector a a =ft> have the form: 

nN d 

M(L, K) = / [] dVidfji f) h r) kl ■ ■ ■ f}i m T] km e^ 1 , (15) 
J i=i 

with L = {li, ■ ■ ■ , l m } and K = {ki, ■ ■ ■ , k m }. The products fjrj are ordered in such a way 
that h < h < ■ ■ ■ < l m an d k\ < k% < ■ ■ ■ < k m . From reference the result of this type 
of integrals is equal to: 



M(L, K) = (_l)('x+I a +-+U) + (fcl+fc2+-4-fcm) A ( Lj R ^ ( 16 ) 
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where A(L, K) is the determinant of the matrix obtained from matrix A** by deleting 
the lines {h, ■ ■ ■ ,l m } and the columns {ki, ■ ■ • , k m }. M(L,K) is a co-factor of matrix 
A^ . The Grassmann integrals to be calculated in sector jj are the same type as eq.flI5|). 
Evaluating determinants of non-diagonal matrices of dimension nN d x nN d is still a 
hard task, even if we have restricted ourselves to multivariable integrals of a fixed sector 
oo. Calculating such determinants is a suitable task for computers, and it obviously 
depends on hardware and software resources. Fixing n, for instance, there is an upper 
practical limit for N, so that the calculation of determinants is feasible. One possibility 
for evaluating eq.(14) is that of assigning different values for N and, from the results 
obtained, trying to extrapolate for an arbitrary value of N. If we are lucky, some recursion 
expression for eq. fll4|) for all N, could be recognized. 

Our approach to calculate the integral (|T5|) is, for fixed n and arbritary N, to explore 
the block-structure of matrices A " ", o —\ and o —I, diagonalizing it through a similarity 
transformation 



where the matrix D is, 

^N d xN d ®N d xN d 
_ ®N d xN d ^2^-N d xN d 

V ®N d xN d ®N d xN d 

Aj, i = 1,2, •••,7i, are the eigenvalues of matrices A " ", o =t, I? and calculate the co- 
factors of the matrix D. The j^-column of matrix P is the eigenvector of A aa associated 
to the eigenvalue Xj. Each eigenvalue of matrix A aa has degeneracy N d . The matrices are 
not hermitian, thus some eigenvalues are complex. In the following, we will be working 
on the oo =ff sector; however, the results for the oo =H sector are analogous since 
A TT = A u . 

We will apply the following transformation of variables, 



®N d xN d \ 
®N d xN d 

X n ^N d xN d J 



7] ' = p- 1 r ] and rf = r/P, (19) 

where rj 1 = {rj[, • • • , i]' nNd } and ff = {fj 1 ^ ■ ■ ■ , ff nNd \. The jacobian of the transformation 
(|19l) is equal to one. 

Due to the fact that A^ is a block-matrix, the matrix P also has 
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a block structure. This fact implies that transformations ( p!9|) do not mix up lattice 
indices. 

In a schematic way, the integrals M(L,K) (eq.flI5|)) become: 



- nN d 

M(L,K) = / UdTHdfk (fjP- 1 ) h (Pr 1 ) kl ---(fjP- l ) lm (Pr 1 ) 

8=1 



J2 ViDijVj 



(20) 



where -Djj are the entries of the diagonal matrix D. The expression M(L, K) fits into the 
form of eq.flTB]), and hence corresponds to some co-factor of the diagonalized matrix D 
(eq.( P^) ). It is very simple to calculate these co-factors, and the matrix P is the same for 
any self-interacting fermionic model. 



In Appendix |A| we present the derivation of the eigenvalues and eigenvectors of matrix 
D for arbitrary values of n and N. From eqs. (|65|) and fl66|), for arbitrary value of n, we 
have that 



and 



with v. v' = 0, 1, 



and 



V, 



X Z7T (Q,,f 

e~ 



n 



» 



1 



n 



, n — 1 , and 

/ ( ra )n 

' POO 1 N d xN d 



\ ( n ) 11 

\Pn-l,Q JL N d xN d 



e _ffi(2i/+l)(i/+l) 



.-1 



N d xN d 



\ ( n ) 11 

\ 9n-l,0 Jl ./V< i xAr <i 



( n ) 11 \ 

Pa,n-l 1 N d xN d 



( n ) n / 

Pn-l,n-l A N d xN d I 



( n ) n \ 
g ,n-l 1 iV £i xA rd ' 



11 / 

yn-Kn-l^A^xiV 1 * / 



(21) 
(22) 



(23) 



(24) 



The diagonal elements of matrix D are: 



^ = 0,1, 



,71 — 1. 



(25) 



where the eigenvalues are iV d -fold degenerated, N d being the number of lattice sites. Due 
to lattice translation symmetry, we should note that the elements p^J, and q„J, do not 
carry any lattice site index. 
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This is a general approach, and it can be applied to any self-interacting fermionic 
model. The important point here is that the relations (|2l|)-(|2"5|) are valid for any self- 
interacting fermionic model with space translation symmetry. 



3 Application to Hatsugai-Kohmoto Model 

The calculation of an exactly soluble model is a nice way to test a new approach. Hatsugai 
and Kohmoto ||To| proposed a toy model (HK model) that shares the atomic and band lim- 
its of the Hubbard model |TT|. Using the Green function and path integral approaches, 



Nogueira and Anda|l2| established the equivalence of this model (with unrestricted hop- 
ping) and the Hubbard model (with infinite-range hopping). 

In this section we derive the grand canonical partition function of the HK model using 
the results presented in section 2. The hamiltonian of the HK model in momentum space 
is PI 



H = EE £(fcK(£) + f/E n #)n#) 

k °"=1U k 

= EH(4 (26) 

k 

where n a (k) = aJ a (k)a. a (k) and aj.(/c) (a (T (fc)) is the creation (destruction) operator of 
an electron with momentum k and spin o. The function e{k) = —2tJ2i=i cos ki, k = 
(ki, k 2 , fc 3 ), corresponds to the nearest hopping of the electrons in the dual-space lattice. 
U is the strength of the repulsion between electrons with the same momentum k but 
opposite spin conponents. 

From eqs. (|l|) and (26), the grand canonical partition function of the HK model is 



m?) = Tr[l[e-^} 

k 

= U[Tr,e-^]. (27) 



We have 



K(k) = E n a (k) + Un^n^k), {21 

<r=hl 
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where we define: A(k) = e(k) — fi, and \x is the chemical potencial. In eg. (p7|) , the symbol 
Tr % stands for the trace for a fixed vector k, whereas Tr represents the trace for all fc's. 
The high temperature expansion for the grand canonical partition function Z(j3; fi) is 

Tr u [e-^] = ^i//?«Tr[K"(fc)]. (29) 

n=0 U - 

Since all the operators on the r.h.s. of eq.(|28|) commute, we can apply the Newton's 
multinomial expression to write Tr.[K n (/c)] as 



— * ' _ 77 I —t — » — » 

Tr,[K n (k)}= J2' j i A" 1+ " 2 (fc)t/" 3 Tr e [n^ 1+n3 (A;)nf +n3 (A;)]. (30) 

ni,n 2 ,n 3 =0 n l- n 2- n 3- 

The symbol X)' means that the summation indices satisfy the condition rii + ri2 + = n. 

Let lx, I2 and Z3 be the integers that determine the lattice vector k. The mapping (|13|) 
takes the index k into the index L = 1% + — 1)N + (/3 — l)iV 2 , where N is the number 



of points in the momentum lattice in each direction. In the sum on the r.h.s. of eq. (|30| ) . 
we calculate the trace for a fixed k, which means that in eq.(6) we take a single point in 
the momentum lattice (N =1). Then, 

Tr u [n^(k) nf+-(£)] = 1^ x 1^ (31) 

where 

n 

r « E 9x(£)A^,,j(i) 

4I.«3 = / II d-ni(L)dfii(L) e'^ x 
J i=i 

xfj (L)r]o(L) ■ • •77 ni _ 1 (L)?7 m _ 1 (L) fj ni+n2 (L)r) ni+n2 (L) ■ • • r/ n _ 1 (L)?7 n _ 1 (L), (32) 



and 



2n 

/• 2 « E %(i)A| i / r?j(L) 

J nin 3 = / II dr)j(L)dfjj(L) e'^ x 

J J=n+1 

(L)x 

x fj n+ni+n2 (L)r] n+ni+n2 (L) ■ ■ ■ f]2n-i(L)r]2n-i(L) . (33) 

The matrices A " ", cr =|, I, are given by eq.(ll) with N = 1. According to eqs. fll5l) 
and (|T6|) , the presence of rfs (and r/'s) in the integrand on the r.h.s. of eqs. (^) and 
(j33[) allows one to evaluate the integrals as the determinants of matrices obtained after 
deletion of lines (and columns) of the matrices A " ", a =f, j. For this particular model, it 



10 



turns out easier to apply eq. (|16|) directly, rather than using the similarity transformation 
(|T?|), since the lattice is unidimensional. We should mention that for iV = 1, we recover 
the case of the anharmonic fermionic oscillator, which has been considered in a previous 
work 0. Now we discuss the values of X^j n3 (eq. fl32|)), in view of the possible values of 
the indices (rii, n 2 , n 3 ). 

i) rii = n, n 2 = and n 3 = 0. 

In this case the first n lines and the first n columns of matrix are deleted; hence, 



4!o = 1- (34) 



it) n\ — 0, n 2 = n and n 3 = 0. 

In this case no lines or columns are deleted in A^; so, 



J TT = det(A TT ) = 2. (35) 



n. 



in) ni = 0, ri2 = and n 3 

This case is equal to case % and therefore 



Hi = !■ (36) 

iv) rii ^ 0, n 2 7^ and n 3 ^ 0. 

In this case, the first n% lines and columns are deleted, as well as the last n 3 lines and 
columns of matrix A". The triangular matrix thus obtained has its determinant equal 
to 1, for any value of n. Then, 

ZlU = 1- (37) 

Equivalent results are valid for X** ns . 

From the results (|34l)-([37D and the equivalent results for n , we have 

Tr^^ik) nf +n \k)\ = (1 + 5„ 1+n3i o) (1 + 5„ 2+ „ 3 ,o), (38) 
that substituted in eq. (|30|) gives 

Tr t [K n (k)] = [2A{k) + U] n + 2A n {k). (39) 
Returning to eqs. (p?D and (]29|) , we finally get 
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Z(J3; fj.) = nil + e-^^ + 2e~ p ^% (40) 

k 

that gives the same free energy density found in reference \Y1\ . 

4 Conclusions 

Calculations involving fermionic fields do demand some extra care, in comparison to the 
manipulation of bosonic fields. For this reason, fermionic models are usually bosonized, 
in a strategy designed to avoid the "annoying" fermionic features. However, moments 
of grassmannian multivariable integrals can be easily calculated, as shown in eqs. (|15D 
and (0). In this paper we have presented a new approach, based on the explicit use of 
Grassmann algebra properties, to the problem of calculating the coefficients of the high 
temperature expansion of the grand canonical partition function for any rf-dimensional 
self-interacting fermionic model (d = 1,2,3, •• ■). We have explored the results (|T5D and 
(p"6|) and the possibility of performing the similarity transformation (|T7|) for a system 
with arbitrary dimension d and arbitrary number of lattice points N d . It is important to 
point out that the matrices A aa (cr =t, I) are model-independent; they are solely related 
to kinetical aspects of the approach. To simplify the notation, we considered that the 
number of points in each direction of the lattice is the same, but the results derived are 
still valid if this is not true. The fact that our results are analytical allows us to obtain 
the thermodynamical limit for any self-interacting fermionic model. 

As a simple example of application of the method (that does not explore all of its 
features, though), we have considered the Hatsugai-Kohmoto model, which is diagonal in 
momentum space and had been solved by other approaches. We have derived its grand 
canonical partition function, and obtained the same free energy density found in the 
literature [|T5 . 



The most important features of this method appear when the hamiltonian has non- 
commuting terms and, consequently, Newton's multinomial expansion does not apply. 
Equations (|ll^-(|T6l) and the similarity transformation (|17]) then become the keystone of 



our analytical results. That is the case of the Hubbard model fLl| ; the grand canoni- 
cal partition function for the unidimensional version of this model is known in integral 
form|13[. A closed expression for this function was obtained by Takahashi fll3 |, within 



certain limits only. By the application of the approach we have presented here, we are 
currently calculating the coefficients of that partition function, up to order /3 5 and for 
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any value of the parameters of the model, as well for any value of the chemical potential. 
These calculations will soon be submitted to publication. 
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Appendix 



A Calculation of Eigenvalues and Eigenvectors of Ma- 
trix A aa 

This appendix is devoted to calculating the eigenvalues and eigenvectors of the matrix 
A acr , defined in eq.(|ll|), as well as determining the matrices P and P -1 that diagonalize 
it (see eq.(ff7I)). 



The characteristic equation for A arj is 
/ (1 — A) l N d xN d —t N d xN d 



det 



d xN d 



®N d xN d 

t N d xN d 
= 0. 



1 — A) t N d xN d 
Q) N d xN d 

®N d xN d 



Qv d xN d 
— l N d xN d 

Q) N d xN d 
®N d xN d 



Q) N d xN d 
Q) N d xN d 

T — A) l N d xN d 
Q) N d xN d 



Q) N d xN d 
Q) N d xN d 



-1 



N d xN d 



1 — A) t N d xN d J 
(41) 



Observe that this matrix (of total dimension nN d x nN d ), consists of a n x n block matrix, 
each block having dimension N d x N d . Moreover, these blocks are either null matrices 
Q N d xN d or proportional to the identity matrix t N d xN d. 

We will demonstrate a useful property of the determinant of a block matrix in which 
all blocks are diagonal, such as the previous matrix. Take a block-matrix M composed of 
blocks B^l, namely 
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M 



/EM B^ 

B [2 ' 1] B [2 ' 2] 



(42) 



VbI 71 ' 1 ' B' n ' 2 ] ■ ■ ■ B[ n,n l / 



where each block B^ is diagonal: 



(43) 



where a, /3 = 1,2, ...,N d . (No summation over repeated indices is implied). We define a 
"determinant-like" matrix function F upon the blocks B^ as 

F = £ 9tt BWBM...BW, (44) 

e 1 ,e 2 ,...,6»„=i 

where ee^a,. e n ^ s the Levi-Civita symbol in n-dimension. Obviously, F and the blocks 
B^'l have all the same dimensions, N d x N d . Using eq . (f43|) , we have 

n 

F ab = 5 a , £ se lA ,...,eM' 6l] ^ M ■ ■ ■ ^ M )- (45) 

6>i,0 2 ,-,0n=l 

From eq.(^5|) and the definition of the determinant, we obtain 



det F = det T x det T 2 . . . det r n , 



(46) 



where we have defined n matrices T p of dimension N x N as 



T 1 = b 1 "''" 1 

V 1 - pj™ — "p 



so that 



(47) 



N 

det r p = ^,..,«« H ,UJl] h v M ■ ■ ■ bjr lM, - i] H n, "" ] ) (48) 

Thus, the evaluation of det F is equivalent to the evaluation of the determinant of a block 
matrix T, defined as 



xN d 

®N d xN d r 2 



'N d xN d 



xN d 



®N d xN d \ 

Q) N d xN d 



(49) 
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However, T and M only differ by an even number of permutations of lines and columns! 
More specifically, M can be recovered from T if we reorder the lines of the latter according 
to the pattern 

(1,2,. ..,nN d ) — > (1, N d + 1, 2N d + 1, . . . , (n — l)N d + 1, 

2, N d + 2, 2N d + 2, (n - l)iV d + 2, . . . 

. . . , N d — 1 , 2iV d - 1, 3N d - 2, . . . , niV d - 1, 

N d , 2N d , 3N d ,..., nN d ) (50) 

- i.e., the 1 st line is left untouched, the 2 nd line is replaced by the (N d + l) th line, etc.- and 
then have the columns of the resulting matrix reordered in the same fashion. (The same 
result is obtained if we reorder columns before lines.) As the total number of permutations 
is even, we have 

detM = detr. (51) 
Combining fl5lD , (f46|) and ((44]), we finally obtain 

n 

det M = det F = det( V e 6l ,e 2 ,...,e n . . . B^) (52) 

e u e 2 ,...,e n =i 

In conclusion, if the matrix M is composed of diagonal blocks B^l, the determinant of 
M is equal to the determinant of the matrix F, defined as a "determinant-like" function 
upon the blocks B^l 



Turning our attention back to eq.(^Tj), we expand the determinant in terms of "cofac- 



tors" , based on the last "line" of blocks: 



det (-1) 1+ ™ l NdxNd (-ltfd^d)"- 1 + {-l) n+n (1 - A) t NdxNd (1 - A)"- 1 J = 

(53) 

which yields the characteristic equation 

(1 + (1 - \) n ) Nd = 0. (54) 
There are n distinct eigenvalues A&, each one with multiplicity N d , given by 

A fc = l-e^+D, (55) 
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where k — 0, 1, 2, ...,n— 1. Observe that if A& is an eigenvalue, so is its complex conjugate: 
A£ = A n _fc_i. Let us denote by an eigenvector of A aa associated to A&. It has the 
structure 



Vv n / 



(56) 



where each Vj, % = 1, 2, n is a 1 x N d matrix. We obtain 



Vj = — (1 — \ k y £, where % — 1, . . . , n — 1 (57) 

and £ = v n is an arbitrary column vector of dimension N d . There are N d possible linearly 
independent choices for £, corresponding to N d distinct eigenvectors associated to the 
same eigenvalue A fc . We choose them to be 



Sfc — 






V o / 

r(0 



t(2) 



/ \ 

-1 



V o J 



AN d -l) 



/ \ 


-1 

V o J 



{N d ) 



( o \ 




v-iy 



(58) 



so that each ^} corresponds to an eigenvector Vj^, where / = 1,2, ...,N d , associated to 
the eigenvalue A&. 

The matrix P that diagonalizes A aa can be obtained by concatenating all eigenvectors 
VjjP for all eigenvalues A&, k — 0, 1, n — 1, up to a normalizing factor R: 



^ Pbo^-N d xN d 



( n ) -n 

\Pn-l,0 R N d xN d 



( n ) 11 \ 

Po,n-l R N d xN d 

( n ) H 
Pn-Kn-l-WxA^ / 



(59) 



where 



with v, i/ — 0,1, 



V 



(n) 



) =jRe f(2,' + l)(. + D ) 



■ , n — 1 . The matrix 
/ ( n ) n 

» <?00 JLjVxiV 



P 1 



(n) -n 
\ 9n-l,0 J WxJV 



( n ) 11 \ 

%,n-l*NxN 

(n) -|, 

Qn—l,n—lNxN / 



(60) 



(61) 



16 



where 



q^l = B! e~^' +1 ^ +l \ (62) 
is the inverse of P, upon a suitable choice of R and R'; i.e, 

R — R' — (63) 

so that they satisfy the relation 

n-l ^ ^ 



P=0 



Hence, 



p W = JLeSW+iXr+Q (65) 



g H = J_ e -^(2,' + i)(, + i)_ (66) 
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